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ABSTRACT 

We report on optical high-dispersion spectroscopic monitoring observations of the 
Be/X-ray binary A0535+26/V725 Tau, carried out from November 2005 to March 
2009. The main aim of these monitoring observations is to study spectral variabilities 
in the Be disc, on both the short (a week or so) and long (more than hundreds of days) 
timescales, by taking long-term frequent observations. Our four-year spectroscopic 
observations indicate that the V/R ratio, i.e., the relative intensity of the violet (V) 
peak to the red (R) one, of the double-peaked Ha line profile varies with a period 
of 500 days. The H/3 line profile also varies in phase with the Ha profile. With these 
observations covering two full cycles of the V/R variability, we reconstruct the 2-D 
structure of the Be disc by applying the Doppler tomography method to the Ha and H/3 
emission line profiles, using a rigidly rotating frame with the V/R variability period. 
The resulting disc structure reveals non-axisymmetric features, which can be explained 
by a one-armed perturbation in the Be disc. It is the first time that an eccentric disc 
structure is directly detected by using a method other than the interferometric one. 
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1 INTRODUCTION 

Classical Be stars (Be stars for short) are non-supergiant B- 
type stars which show or have shown Balmer lines in emis- 
sion. The striking feature of the Be stars is their high rota- 
tion speeds. They rotate at a speed close to critical so that 
the surface gravity mostly balances with their centrifugal 
force around the equator. The mass ejected from the pho- 
tosphere forms an circumstellar envelope, called a Be disc, 
which is another striking feature of Be stars. The Be stars 
exhibit complicated line profiles containing an absorption 
component from the photosphere and an emission compo- 
nent from the disc. Other observed features of the Be stars 
are IR-excess due to /-/ and f-b emission, and linear polar- 
ization due to electron scattering in the Be disc. The emis- 
sion line profiles, reflecting the state of the Be disc, show 
many kinds of variability on the timescales from days to 
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several decades, e.g., disappearance and reformation of the 
disc, and line profile variability (LPV) including V/R vari- 
ability, the variability in the ratio of the intensity of violet 
(V) and red (R) peaks of double-peaked emission line pro- 
files. 

Recently, IPorter fc Riviniusl (|2003l ) reviewed the latest 
results of Be stars: mainly their properties, variabilities, and 
probable mechanisms of the disc formation. The mechanism 
of Be-disc formation is still under debate, because it is diffi- 
cult to model the complex mass ejection process, taking into 
account the various physics at work, e.g., the effect of rapid 
rotation, the radiative force, and non-LTE calculations. The 
common understanding is, however, that the key to Be-disc 
formation is the hi gh rotation speed o f the B e stars [for fur- 
ther discussion, see iPorter fc Riviniusl (|2003h and references 



therein]. Some authors have suggested t hat the high rota- 
tion speed is linked to ste llar evolution (IMcSwain fc Giesl 
120051 : lEkstrom et al~ll2008l ) . 

Be/X-ray binaries are systems which consist of a Be star 
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and a compact object (mostly a neutron star) a nd account 
for a large fraction of high-mass X-ray binaries (|Coe et al. I 
2009). The system has two discs: a Be disc and an accretion 
disc around the neutron star, which is formed by material 
transferred from the Be disc. Hence, the Be/X -ray binaries 
exhibit accretion phenomena as well as Be phenomena. 

Generally, the eccentricity of the Be/X-ray binary is 
not small (> 0.3), which implies that the interaction be- 
tween the two stars and the mass transfer from the Be 
disc to the neutron star de pends on the orbital phase 
ijQkazaki fc Neeueruelal l2001i ). Most of the Be/X-ray bi- 
naries are thus transient X-ray sources. The X-ray out- 
bursts from these systems are divided into two types ac- 
cording to the X-ray luminosity: normal outbursts (type I 
outbursts; ~ 10 36 - 37 erg s~ x at 1 - 20 keV), and giant out- 
bursts (type II outbursts; 10 37 erg s" 1 at 1-20 keV). 
The normal outbursts, which last several days, are known 
to occur around periastron passage of the neutron star. 
lOkazaki fc Neeueruelal (|200lh appli ed the resonantly trun- 
cated Be-disc model developed by iNegueruela fc Okazakil 
|200ll ) to several Be/X -ray binary systems. They concluded 
that the normal outbursts occur in systems with interme- 
diate to high eccentricities where the mass transfer from 
the Be disc to the neutron star takes place at every perias- 
tron passage. On the other hand, the giant outbursts, lasting 
several tens of days with no orbital modulation, is not well 
understood. 

A0535+26/V725 Tau, which was discovered by the Ariel 
V satellite durin g a giant outburst (jRosenberg et al.lll975l : 
ICoe et al. \\w7fb , is one of the best studied Be/X-ray bina- 
ries. The sys tem consists of an X-ra y pulsar of 103 seconds 
spin period dCaballero et al. 1 120071 ) and an 09.7IIIe star 
jGiangrande et al. Ill980h . The recurrence time of the nor- 
mal outbursts and spin-do wn rate of the pulsa r provide the 
orbital period of 111 days dMotch et al. MlQQlh and the ec- 
centricity of 0.47 (|Finger et al. Ill994h . Neither optical pho- 
tometry nor spectroscopy, however, ha s yielded an orbital 
period consistent wit h X-ray data fe.g. JWang fc Gieslll998l ; 
lLarionov et al. Il200ll ). because the optical data shows very 
complex variabilities of line profiles due to the variability of 
the Be star. 

Long-term variabiliti es in this system h av e bee n re- 
ported by many authors. IClark et al. I |l998al lbl. I1999T ) per- 
formed UV, optical and IR spectroscopy and U, B, V- 
band photometry. Seven-year optical spectroscopy exhibited 
quasi-periodic variability with a period of ~ 1 year, and 
fifteen-year photo metry showed variabi lity in phase with the 
spectroscopic one. Hai gh et al. I {2004>) reported cyclic be- 
haviour due to resonant truncation of the Be disc by carry- 
ing out optical spectroscopy and I, J, Tf-band photometry 
for fourteen years. The dramat ic Be disc growth betwee n 
1998 and 2000 was observed by iGrundstrom et al. I (|2007h . 
which also reported V/R variabilities. 

However, no short-term variabilities (a few week or less) 
can be discussed from these observations, since they have 
not been carried out densely enough. Besides, many optical 
spectroscopic observations reported previously had low or 
medium dispersion, and it is therefore difficult to discuss 
small amplitude (< 10 km/s) variability in the disc. 

In this two-paper series, we report on, both of 
the short-term and long-term variabilities of the Be 
disc in A0535+26/V725 Tau system, based on optical 



high-dispersion spectroscopic monitoring observations from 
November 2005 to March 2009. In this paper (paper I), 
we focused on long-term (> years) variabilities. Short-term 
(from days to weeks) variabilities will be discussed in paper 
II (Moritani et al. in preparation). 

In Sec. 2, we summarise our observations. The results 
are given in Sec. 3, and we discuss the long-term variability 
of A0535+26 in Sec. 4. We present our conclusions in Sec. 
5. 



2 OBSERVATION 
2.1 configuration 

We carried out high dispersion optical spectroscopic mon- 
itoring observations of A0535+26 from November 2005 to 
March 2009, mainly at the Okayama Astrophysical Observa- 
tory (OAO) with a 188 cm telescope equipped with HIDES 
(High Dispersion Echelle Spectrograph). Observations were 
also performed at Gunma Astronomical Observatory (GAO) 
with a 1.5 m telescope equipped with GAOES (Gunma As- 
tronomical Observatory Echelle Spectrograph). HIDES cov- 
ers a 1200 A (5500 - 6700 A and 4400 - 5600 A) wavelength 
range until 2007, and 3500 A (4000 - 7500 A) from 2008 
thanks to adoption of mosaicked three EEV 42 - 80 CCDs 
(2048 x 4098 pix 2 ). The wavelength coverage of GAOES is 
1900 A (4800 - 6700 A) and the detector is a EEV 44 - 82 
CCD. 

The typical wavelength resolution R and the signal to 
noise ratio S/N of our OAO /HIDES data around Ha are R 
~ 60000 and S/N ~ 120, respectively. For H/3, R ~ 60000 
and S/N ~ 100. On the other hand, our Ha data obtained 
with GAO/GAOES has R ~ 30000 and S/N ~ 120, and 
the H/3 data has R ~ 30000 and S/N < 100. The exposure 
time for the Ha data with OAO/HIDES ranged from 3600 
to 5400 s, and that of H/3 data was from 3600 to 7200 s. 
With GAO/GAOES the exposure time was 3600 - 9600 s. 
We combined 1200- or 1800- second-exposure data into daily 
averaged spectra. 

The obtained data is reduced in the standard way us- 
ing IRAE0 echelle package - subtraction of bias, flat fielding, 
calibration of the wavelength using Th-Ar lines, normaliza- 
tion of the continuum, and correction to the helio-centric 
orbit. 



2.2 observation periods 

Our observation log is listed in table [1] The observation pe- 
riod in column 1 is described below. From columns 5 to 8 
are listed the orbital phases corresponding to three different 
ephemerides by different authors and the X-ray light curve. 
The time of phase (HJDo) is given by the following equa- 
tion, 

HJD() = H J D or i g i n + P or lj x E, (1) 



1 IRAF (Imaging Reduction and Analysis Facility) is a software 
system for the reduction and analysis of astronomical data, sup- 
ported by NOAO. 
|http://iraf.noao.edu/| 
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Table 1. Observation log. Observations are divided into nine periods displayed in column 1 (see text for details). Column 3 shows HJD 
of the mid exposure time of each observation. Orbital phases corresponding to three different ephemerides are listed in column 5 to 
column 7, while the phase based on the X-ray light curve is given in column 8 (for the definition of each ephemeris, see tabled- Columns 
9 and 10 show whether the Ha and H/3 line profiles respectively were observed (Q) or not ( — ). 



Observation 




















Period 


Date 


HJD 


Observatory 


(Ephemeris 1) 


(Ephemeris 2) 


(Ephemeris 3) 


(Outburst) 


Ha 




1 


2005 Nov. 24 


2 453 699.060 


OAO 


0.576 


0.066 


0.773 


0.727 


o 






2005 Nov. 25 


2 453 700.060 


OAO 


0.585 


0.075 


0.782 


0.736 


o 






2005 Nov. 27 


2 453 702.060 


OAO 


0.603 


0.093 


0.801 


0.754 


n 






2005 Nov. 29 


2 453 704.031 


OAO 


0.621 


0.111 


0.818 


0.772 








2005 Nov. 30 


2 453 705.081 


OAO 


0.631 


0.120 


0.828 


0.782 








2005 Dec. 01 


2 453 706.075 


OAO 


0.639 


0.129 


0.837 


0.791 








2005 Dec. 03 


2 453 708.651 


OAO 


0.663 


0.153 


0.860 


0.814 


o 


— 


II 


2006 Dec. 18 


2 454 088.100 


OAO 


0.072 


0.593 


0.310 


0.256 






III 


2007 Nov. 07 


2 454 412.184 


OAO 


0.984 


0.531 


0.256 


0.196 


o 


o 




2007 Nov. 08 


2 454 413.172 


OAO 


0.993 


0.540 


0.265 


0.205 




p 




2007 Nov. 09 


2 454 414.141 


OAO 


0.001 


0.549 


0.274 


0.214 


n 


p 




2007 Nov. 10 


2 454 415.232 


OAO 


0.011 


0.559 


0.284 


0.224 


p 






2007 Nov. 11 


2 454 416.166 


OAO 


0.019 


0.567 


0.292 


0.232 


p 


p 




2007 Nov. 13 


2 454 418.159 


OAO 


0.037 


0.585 


0.311 


0.250 


p 


p 




2007 Nov. 14 


2 454 419.238 


OAO 


0.047 


0.595 


0.320 


0.260 


p 


p 


IV 


2007 Dec. 16 


2 454 451.168 


GAO 


0.334 


0.885 


0.611 


0.550 


o 


o 




2007 Dec. 19 


2 454 454.294 


GAO 


0.362 


0.913 


0.639 


0.578 


p 


p 




2007 Dec. 20 


2 454 455.264 


GAO 


0.371 


0.922 


0.648 


0.587 


p 


p 




2007 Dec. 26 


2 454 461.090 


GAO 


0.423 


0.975 


0.701 


0.640 


p 


p 




2008 Jan. 02 


2 454 468.022 


OAO 


0.485 


0.037 


0.764 


0.702 


O 


o 




2008 Jan. 31 


2 454 497.004 


OAO 


0.746 


0.300 


0.027 


0.965 


o 


o 


v 


2008 Mar. 11 


2 454 537.002 


GAO 


0.105 


0.663 


0.391 


0.328 


p 






2008 Mar. 21 


2 454 547.023 


OAO 


0.195 


0.754 


0.482 


0.419 


p 


p 


VI 


2008 Oct. 01 


2 454 741.208 


OAO 


0.940 


0.514 


0.247 


0.180 


o 


o 




2008 Oct. 11 


2 454 751.162 


OAO 


0.029 


0.604 


0.338 


0.271 


o 


o 




2008 Oct. 11 


2 454 751.229 


GAO 


0.030 


0.605 


0.338 


0.271 


o 


o 




2008 Oct. 15 


2 454 755.186 


GAO 


0.065 


0.641 


0.374 


0.307 


o 


o 




2008 Oct. 19 


2 454 759.141 


GAO 


0.101 


0.677 


0.410 


0.343 


o 


o 




2008 Nov. 04 


2 454 775.102 


GAO 


0.244 


0.821 


0.555 


0.488 


o 


o 




2008 Dec. 12 


2 454 813.068 


GAO 


0.586 


0.166 


0.901 


0.832 


o 


o 


VII 


2008 Dec. 25 


2 454 826.207 


OAO 


0.704 


0.285 


0.020 


0.952 


o 


o 




2008 Dec. 26 


2 454 827.097 


OAO 


0.712 


0.293 


0.028 


0.960 


o 


o 




2008 Dec. 27 


2 454 828.129 


OAO 


0.721 


0.302 


0.038 


0.969 


o 


o 




2008 Dec. 28 


2 454 829.187 


OAO 


0.730 


0.312 


0.047 


0.979 


o 


o 


VIII 


2008 Dec. 29 


2 454 830.165 


OAO 


0.739 


0.321 


0.056 


0.987 


o 


o 




2008 Dec. 30 


2 454 831.157 


OAO 


0.748 


0.330 


0.065 


0.996 


o 


o 




2008 Dec. 31 


2 454 832.112 


OAO 


0.757 


0.338 


0.074 


0.005 


o 


o 




2009 Jan. 01 


2 454 833.008 


OAO 


0.765 


0.346 


0.082 


0.013 


o 


o 




2009 Jan. 02 


2 454 834.021 


OAO 


0.774 


0.356 


0.091 


0.022 


o 


o 




2009 Jan. 03 


2 454 835.028 


OAO 


0.783 


0.365 


0.100 


0.032 


o 


o 




2009 Jan. 05 


2 454 837.114 


OAO 


0.802 


0.384 


0.119 


0.050 


o 


o 




2009 Jan. 06 


2 454 838.138 


OAO 


0.811 


0.393 


0.129 


0.060 


o 


o 




2009 Jan. 07 


2 454 839.180 


OAO 


0.820 


0.402 


0.138 


0.069 


o 


o 




2009 Jan. 08 


2 454 840.194 


OAO 


0.829 


0.412 


0.147 


0.078 


o 


o 




2009 Jan. 10 


2 454 842.186 


OAO 


0.847 


0.430 


0.165 


0.096 


o 


o 




2009 Jan. 12 


2 454 844.157 


OAO 


0.865 


0.447 


0.183 


0.114 


o 


o 


IX 


2009 Mar. 12 


2 454 902.965 


GAO 


0.393 


0.981 


0.718 


0.648 


o 


o 
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Figure 1. Power spectrum of the RXTE/ASM data. The ar- 
row in the spectrum indicates the most plausible value, which is 
corresponding to 110.24 days. The highest peak corresponds to 
366.45 days, introduced by yearly sun angle constraints in RXTE 
observations of A0535+26. 

where H JD or i g i„ is the origin, P or b the orbital period, and 
E an integer indicating the cycle number. H J D or i g i n cor- 
responds to a periastron passage in the case of ephemeris, 
while it corresponds to the time when the X-ray flux starts 
to rise in the case of X-ray light curve (see below). We used 
three different ephemerides reported by different authors be- 
cause of difficulty in determining the time of periastron pas- 
sage. The H J D or igi n and P OT b of each ephemeris are sum- 
marised in table [2] 

The outburst phase listed in column 8 is obtained by us- 
ing RXTE/AS1V0 one-day averaged data from JD 2450133 
(2 February 1996) to JD 2454881 (19 February 2009) as fol- 
lows: First, we determine the orbital period by the Fourier 
analysis of the X-ray light curve. To remove uncertainty due 
to giant outbursts, the data between —1.5 and 4.0 ASM 
Unit counts s^ 1 (1 Crab is approximately 75 ASM Unit 
counts/s at 2 - 10 keV) were used for our period analy- 
sis. The power spectrum is shown in figure [T] The arrow 
in the power spectrum indicates the most plausible peak: 
0.00907±0.00005 [day -1 ], corresponding to P or6 =110.2±0.6 
days. Other peaks are overtones or harmonics with one year. 
Then, the beginning time of an X-ray normal outburst, HJD 
2453398.4, is taken as H JD or i 9 i n . The obtained period and 
the origin are listed in table [2] together with those in the 
three different ephemerides. Figure [5] shows the X-ray light 
curve folded with these P or b and H J D or iain by the n x m 
bin method used in lCoe et al. "l (|2006h . The light curve shows 
that normal outbursts reach the maximum flux at (j>x = 0.05 
- 0.08. 

Our monitoring term is divided into nine periods (as 
shown in column 1 of table [TJ. We carried out the monitor- 
ing observations, mainly for detecting the variability of the 
Be disc around periastron passage, when the Balmer line is 
expected to change due to the tidal interaction between the 
Be disc and the neutron star, but also for examining long- 
term variabilities such as V/R variations. Observations in 
periods I (in November 2005), III (in November 2007) and 

2 http: / / xte.mit.edu /asmlc / ASM . html 
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Figure 2. X-ray light curve folded by assuming P or [,=110.24 days 
and HJDorigin = 2453398.43, using the n X m bin method in 
ICoe et al. I feOOfl) . The maximal duration of the normal outbursts 
is 4> x = 0.05 - 0.08. 

VIII (in January 2009), around periastron passage with re- 
spect to the 4>2, 0i and 4>x (</>3), respectively, are performed 
densely, aiming for studying short-term variability. Obser- 
vations in other periods are carried out in order to monitor 
A0535+26 between normal outbursts as well as to aim at 
long-term vari ability. 

Recently, iLevine fc Remillardl (|2008t ) reported that 
A0535+26 brightened in X-rays on 12 September 2008 (JD 
2454721) at 21 mCrab at 2 - 10 keV, and that the flux then 
rose up to approximately 48 mCrab (3.6 ASM Unit counts 
s _1 ). Observations in the period VI, from 1 October to 4 
November in 2008, was carried out thanks to the ToO (Tar- 
get of Opportunity) observation program at OAO, in order 
to examine the variability of the Be disc at <f)x = 0.2 - 0.5 
after the normal outburst. 

Based on the ephemeris given bv lCoe et all (|2006l ). the 
observation of periods VII and VIII were scheduled; the next 
periastron passage was predicted to be around 31 December 
2008, since the last outburst occurred on 12 September af- 
ter periastron. The activity reports of the RXTE/ASM team 
declared that normal outburst occurred in the week of 2 - 9 
January 2009; the X-ray photon count started to rise on JD 
2454832. Our observations were performed from 25 Decem- 
ber 2008 thorough 12 January 2009, which is divided into 
two periods; pre-outburst (period VII, until 28 December) 
and around outburst (period VIII, from 29 December). An 
additional observation was made at GAO on 12 March 2009 
(in the period IX). 



3 RESULTS 

Figure [3] and [3] show representative spectra of Ha and H/3 
obtained during the monitoring observations; all spectra will 
be shown in paper II. For reasons of clarity, the spectra are 
shown on the same Y-scale and with linear offsets from each 
other. For three and a half years, both Ha and H/3 line pro- 
files have been in emission, which indicates that the Be disc 
of A0535+26 did not disappear. The profiles, however, have 
changed. In period I, the Ha line profile was double-peaked 
with the violet component stronger than the red component 
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Table 2. The orbital period, P or b, and the origin of the phase, H J D or i g i n , in each ephemeris and those derived from the X-ray light 
curve. See text for detail. 



Phase 


Type 




^orb 


HJ D origin 


references 


<t>l 


Ephemeris 


1 


111.38± 0.11 


2 446 734.3±2.6 


Motch et al. (1991) 


<j>2 


Ephemeris 


2 


110.3±0.3 


2 449 059.2±0.6 


Finger et al. (19941 


4>3 


Ephemeris 


3 


110.0±0.5 


2 450 094±1 


Coe et al. (2006) 


<t>X 


Outburst 




110.2±0.6 


2 453 398.4±1.0 


this work 




6545 6550 6555 6560 6565 6570 6575 6580 
wavelength [angstrom] 



Figure 3. Representative Ha line profiles from our monitoring 
observations. The observation date is annotated on each spec- 
trum. There are offsets between each spectrum for clarity, and 
larger space is inserted between spectra in different periods. Spec- 
tra in period VII are not shown since they were similar to that 
on 29 December 2008. The line profile variability is clearly seen. 
Sec text for detail. 




4850 4855 4860 4865 4870 4875 
wavelength [angstrom] 



Figure 4. Representative H/3 line profiles obtained from our mon- 
itoring observations. As in Fig. \3\ the observed date is annotated 
on each spectrum, and a larger space is inserted between spectra 
in different periods. The H/3 line has kept double-peaked, but the 
V/R ratio had variabilities nearly in phase with Ha line profile 
(see Fig. [5] and text). 



(V > R, the top profile in figure [3]). One year later (second 
spectrum from the top), in December 2006, the profile kept 

V > R, but the V/R ratio has decreased from > 1.7 to 1.22. 
Here, the V/R ratio was determined by smoothing the spec- 
tra and measuring the peak of the intensities normalized by 
the continuum. It turned to V < R (third spectrum from 
the top) in period III (in November 2007). In period IV, the 
Ha line changed from a red-enhanced single-peaked profile 
to a single-peaked profile with a nearly flat top in less than 
six weeks (fourth and fifth spectra from the top), followed 
by the double-peaked profile with V > R after two months 
(in March 2008, beneath the two). From October 2008 to 
March 2009, from period VI to IX, the Ha line profile re- 
mained double-peaked, gradually changing from V > R to 

V < R in more than half a year (sixth spectrum from the 
bottom). Besides, in period VIII, the profile showed an obvi- 
ous variability within less than two weeks. This short-term 
variability, which we think is due to the tidal interaction be- 
tween the Be disc and the neutron star, will be discussed in 
detail in paper II. 

As shown in figure [4] the H/3 line profile was always 
double-peaked as far as our observations were carried out, 
but the V/R ratio exhibited variabilities. In period III, the 
red component was stronger than the violet, i.e. V < R. In 



periods IV and V (from December 2007 to March 2008), the 
V/R ratio of H/3 line changed from V < R to V > R within 
a couple of months. Then, the V/R ratio gradually turned 
again to V > R during the following year. The variability of 
the V/R ratio is almost in phase with that of the Ha line 
profile (see table [3] and figure [SJ. 

The obtained EW and V/R ratio of all the Ha and 
H/3 line profiles are listed in table [3] The typical error in the 
value of EW is ~ ± 0.05 A , and that of V/R ratio is ± 0.03. 
In period IV, the V/R ratio of Ha line profiles could not be 
obtained because the line profiles did not have double peak. 
In the long term, both EW (Ha) and EW (H/3) gradually 
increased in amplitude during the last three and a half years, 
which indicates that the Be disc of A0535+26 has become 
more and more active. 



4 DISCUSSION 

The observed Ha and H/3 emission-line profiles showed vari- 
abilities on timescales from weeks to years. In what follows, 
we focus only on long-term variabilities (> year). Below, we 
adopt the orbital ephemeris <f>x based on the X-ray light 
curve. 
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Table 3. List of the EW and V/R ratio of all the He* and H/3 emission line profiles observed from November 2005 to March 2009. Column 
1 and 2 list the observation date and corresponding HJD. The orbital phase based on the X-ray light curve, tfix (column 8 in table[T}, is 
listed again in column 3. In columns 4 and 5, EW and V/R ratio of Ha are listed, while in column 6 and 7 are listed those of H/3. The 
typical errors in EW and V/R ratio are ~ ± 0.05 A, and ± 0.03, respectively. 



Day 




HJD 


<t>x 


rLW (Ha) 


V/K (Ha) 


T"1I T J / TT Q\ 

-hW (Hp) 


V/R (H^i) 


2005 Nov. 


24 


2 453 699.060 


0.727 


-6.64A 


1.73 




- 


2005 Nov. 


25 


2 453 700.060 


0.736 


-6.72 A 


1.81 


— 


— 


2005 INOV. 


2 ( 


A 45o fUz.UoU 


0.754 


—6.73 A 


1.90 






2005 INov. 


29 


2 453 704.031 


0.772 


—6.68 A 


1.74 






onnc T\T „ 

2005 Inov. 


30 


O (po 7nr noi 

2 453 705.081 


0.782 


—6.85 A 


1.80 






2005 Dec. 


01 


o /ico Tru" nTrr 

2 453 706.075 


0.791 


—6.49 A 


1.80 






2005 Dec. 


03 


2 453 (08.051 


0.814 


— 6.45 A 


1.73 






2006 Dec. 


18 


2 454 088.100 


0.256 


-9.84 A 


1.22 


— 


— 


2007 Nov. 


07 


2 454 412.184 


0.196 


-9.74 A 


0.77 


-0.95 A 


0.94 


2007 Nov. 


08 


2 454 413.172 


0.205 


-8.83 A 


0.76 


-0.84 A 


0.96 


2007 Nov. 


09 


O A A A 1 A 1/11 

A 454 414.141 


0.214 


—9.42 A 


0.79 


—0.85 A 


0.93 


2007 Nov. 


10 


o a cr A a i rr ooo 

A 454 415. 262 


0.224 


—9.05 A 


0.78 






2007 Nov. 


11 


A 454 41D.1DD 


0.232 


—9.55 A 


0.78 


— 0.92 A 


0.93 


2007 Nov. 


13 


o A a A /II o icn 
A 454 41o.loy 


0.250 


—9.31 A 


0.75 


—0.97 A 


0.93 


2007 Nov. 


14 


O A a A /1TOQQQ 

A 454 41y.Zoo 


0.260 


— 9.25 A 


0.72 


—0.99 A 


0.91 


2007 Dec. 


16 


2 454 451.168 


0.550 


-10.45 A 


* 


-1.23 A 


0.70 


2007 Dec. 


19 


2 454 454.294 


0.578 


-9.73 A 


* 


-1.28 A 


0.53 


2007 Dec. 


20 


O A K A A K C OC A 

A 454 455. Zo4 


0.587 


—9.01 A 




—0.95 A 


0.58 


2007 Dec. 


26 


o a c a Acz~t rifin 
A 454 4ul.UyU 


0.640 


—9.32 A 


* 


— 1.06 A 


0.64 


2008 Jan. 


02 


z 454 4oo.U2z 


0.702 


— 10.23 A 


* 


— 1.14 A 


0.84 


2008 Jan. 


31 


2 454 497.004 


0.965 


-10.27 A 


** 


-1.32 A 


1.01 


2008 Mar. 11 


O Atr^A t;Q7 OOO 

2 454 53 1 .002 


0.328 


— 10.66 A 


1.15 






2008 Mar. 21 


2 404 04 ( .023 


0.419 


—9.90 A 


1.36 


—0.93 


1.40 


2008 Oct. 


01 


2 454 741.208 


0.180 


-11.34 A 


1.61 


-1.22 A 


1.73 


2008 Oct. 


11 


2 454 751.162 


0.271 


-10.20 A 


1.56 


-1.40 A 


1.61 


2008 Oct. 


11 


2 454 751.229 


0.271 


-10.51 A 


1.54 


-1.27 A 


1.58 


2008 Oct. 


15 


2 454 755.186 


0.307 


-10.92 A 


1.39 


-1.08 A 


1.42 


2008 Oct. 


19 


2 454 759.141 


0.343 


-11.24 A 


1.30 


-1.08 A 


1.41 


2008 Nov. 


04 


2 454 775.102 


0.488 


-11.83 A 


1.37 


-1.26 A 


1.76 


2008 Dec. 


12 


2 454 813.068 


0.832 


-11.01 A 


1.35 


-1.03 A 


1.34 


2008 Dec. 25 


2 454 826.207 


0.952 


-11.82 A 


1.27 


-1.38 A 


1.16 


2008 Dec. 


26 


9 4^4 897 0Q7 


n Qfin 
u.yuu 


Y2 Q2 A 


1 28 


1 32 A 


1 14 


2008 Dec. 


27 


2 454 828.129 


0.969 


-12.39 A 


1.23 


-1.36 A 


1.28 


2008 Dec. 28 


2 454 829.187 


0.979 


-12.14 A 


1.22 


-1.45 A 


1.07 


2008 Dec. 


29 


2 454 830.165 


0.987 


-12.48 A 


1.20 


-1.39 A 


1.03 


2008 Dec. 


30 


2 454 831.157 


0.996 


-12.44 A 


1.18 


-1.24 A 


1.06 


2008 Dec. 


31 


2 454 832.112 


0.005 


-12.21 A 


1.21 


-1.42 A 


1.15 


2009 Jan. 


01 


2 454 833.008 


0.013 


-12.12 A 


1.22 


-1.35 A 


1.26 


2009 Jan. 


02 


2 454 834.021 


0.022 


-11.74 A 


1.22 


-1.28 A 


1.09 


2009 Jan. 


03 


2 454 835.028 


0.032 


-11.76 A 


1.23 


-1.47 A 


1.08 


2009 Jan. 


05 


2 454 837.114 


0.050 


-12.86 A 


1.10 


-1.41 A 


1.07 


2009 Jan. 


06 


2 454 838.138 


0.060 


-11.67 A 


1.07 


-1.32 A 


0.98 


2009 Jan. 


07 


2 454 839.180 


0.069 


-12.46 A 


1.05 


-1.54 A 


1.01 


2009 Jan. 


08 


2 454 840.194 


0.078 


-11.45 A 


1.05 


-1.25 A 


0.96 


2009 Jan. 


10 


2 454 842.186 


0.096 


-12.19 A 


1.02 


-1.32 A 


0.88 


2009 Jan. 


12 


2 454 844.157 


0.114 


-12.91 A 


1.00 


-1.47 A 


0.78 


2009 Mar. 12 


2 454 902.965 


0.648 


-12.61 A 


0.84 


-1.74 A 


0.65 



* : Red-enhanced profile 
** : Single-peaked profile 
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JD - 2400000 oscillation phase fl> 



Figure 5. V/R ratio of the Hq (filled circles) and H/3 (trian- 
gles) line profiles vs. HJD in our observation in periods I - IX. 
A remarkable similarit y is seen in the variation p atterns of both 
lines. The V/R ratio in lGrundstrom et al. I ||2007T ) is also plotted 
(squares). The solid line is a sine curve with the period of 500 
days (see also Fig.0. 






Ml V ^11 

\l iJ v , v V v * \i , v i v u 

" 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 
frequency [day' ] 

Figure 6. Power spectrum of the V/R ratio using the data from 
HJD 2453500 to 2455000. The arrow in the spectrum indicates 
the most plausible value, which corresponds to 500 days. 



4.1 The V/R Variation 

As mentioned in the previous section, the relative strength of 
the violet component of Ha compared to the red component 
was greater than 1 in November 2005 and in December 2006. 
Figure 2 of iGrundstrom et al. I l|2007t ). however, shows that 
V/R < 1 in October 2006 (square marks in Fig. [5]). These 
results suggest that the V/R ratio turned from > 1 to < 1 
and then back to > 1 again in 2006 (around JD 2454020). 
The V/R ratio also turned from < 1 to > 1 in February 
2008 between periods IV and V (around JD 2454500). It 
then became < 1 again between periods VIII and IX (around 
JD 2454850). This implies a (quasi-)periodic variation of the 
V/R ratio, which may be induced by a o ne-armed oscillation 
|Okazakilll99ll ; iPapaloizou et al. II 19921 ) . 

We examine the long-term variation of the V/R ratio 
and EW of Ha line from 1979 to 2009 and that of H/3 from 
1975 to 2009, using the data fr om this work, with thos e 
available from previous studies [Hutchings et al. I l| 19781 ) . 



Figu re 7. V/R. variation of Ha since 2005 l|Grundstrom et al. I 
2007), and our observations). The data are folded on 500-day pe- 
riodicity. The phase is defined as the time when the V/R ratio 
turns from > 1 to < 1. The solid line is a simple sine curve with 
the same period. 



Aabl dl985l), IClark et al. I Jl998al) , lHaigh et al. I l|2004h . 
Coe et al. lipOOrj) . and lGrundstrom et al. I (120071 )]. The V/R 



ratio of the Ha seems to have been varying in the range be- 
tween 0.7 to 2, while that of Kf3 between 0.5 and 1.9. The 
V/R ratios of both lines are almost in phase, although these 
lines were seldom observed simultaneously in the previous 
studies. We apply the Fourier analysis for the determina- 
tion of the period of the V/R variation. Here, we restrict 
our analysis to the data of Ha lin e profiles obtained between 
HJD 2453500 to HJD 2455000, bv lGrundstrom et al. I (|2007D 
and this work, because the data are dense enough to derive 
a reliable period. Observations before the 2005 giant out- 
burst, which incl ude dense observations from HJD 2451000 
to HJD 2452000 ijGrundstrom et al. 1120071 ; ICoe et al. Il2006l ; 
IClark et al.~l ll998a), are excluded from the analysis, since it 
is quite possible that the Be disc have reformed after the 
outburst. 

The period of the V/R variation is determined to be 
500±15 days (figure [6| at the 90 % confidence level via 
X 2 test, which is sig nificantly longer than that reported in 
IClark et al.l l| 1998a) ) (~ 1 year). The V/R ratios of the Ha 
spectra from HJD -2453500 to HJD -2455000 folded on 
500 days of period are plotted in figure [7J where the ori- 
gin of the V/R variation phase &v/r is taken at the epoch 
when the V/R ratio turns from > 1 to < 1. We note that 
a cyclic behaviour of the V/R ratio is clearly seen with this 
period. The duration when V > R is longer than that for V 
< R. The simple sine curve wit h this period is als o drawn 
in figure [7] as well as in figure 151 iReig et al~1 (|2005l ) studied 
the period of V/R variations for several Be/X-ray binaries 
(in their table 3). For A0535+26/V725 Tau, the V/R pe- 
riod was suggested to be 1 - 1.5 years l|Haigh et al. 1 120041 ; 
IClark et alll998al ). which is consistent with ours. We need 
further observations for a more accurate period, especially 
bet ween <b v / R = 0-4 and $ v /r = 0.7 (see figure 0. 

lOktariani fc Okazakil (2009) studied one-armed oscilla- 
tions of Be disc in circular bin aries, tak i ng in to account a 
three-dimensional effect, which lOgilviel (2008) found pro- 
vides an important contribution to the confinement of the 
oscillations to the inner part of the Be disc. They calculated 
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that the period of the fundamental mode of the one-armed 
oscillation in a BOV-type Be star ranges from 1.4 years to 
2.3 years, which depends on the binary separation D. Given 
that the mass Mi and radius Ri o f the Be star V725 Tau ar e 
respectively 14 Mq and 15 Rq (|Giangrande et al~lll980l ). 
and that the mass M2 is typical as a neutron star (~ 1 
Mq), D/Ri is approximately 17 in A0535+26, where D is 
estimated by the Kepler's third law, 



3 _ G(Mi + M 2 ) p2 

U - J~2 ^orb- 



(2) 



If we neglect that the luminosity class of the Be star in 
A0535+26 is different than in the i r calcu lations, we see from 
figure 3 of lOktariani fc O kazaki (2009) that the period of 
the one-armed oscillation in this system is ~ 2 years, which 
roughly agrees with our result. However, the period of a one- 
armed oscillation is, in general, sensitive to small changes in 
stellar and disc parameters, and estimated uncertainty of 
the oscillation period is approximately a factor of two. More 
sophisticated calculations, therefore, are needed for further 
comparison with the observation. 

The variation of EW(Ha) and EW(H/3) is given in figure 
[8] The Ha line has been in emission for m ost of the last thirty 
years ; only in 1994, it was in absorption (|Grundstrom et al. I 
120071 ) . The growth of the Be disc is correlated with EW. The 
larger values of jEWj correspond to denser and/or larger Be 
discs. Regarding the las t ten years, a giant ou tburst occurred 
in 2005 (JD=2453507, iTueller at al. 1 120051 s ) . around which 
the emission from the Be disc decreased. Afterwards, the Be 
disc has kept growing. 



4.2 The Structure of the Be Disc 

In order to investigate the two-dimensional structure of the 
Be disc in A0535+26, we carried out the Doppler tomogra- 
phy method using Ha and H/3 spectra. The Doppler tomog- 
raphy method translates LPV into a distribution map of the 
emission component. That is, LPV along the periodic phase 
(usually the orbital phase) can be projected on a 2-D veloc- 
ity plane (V x , Vy) corresponding to the Doppler shift of the 
line profile. This method is often used for circular binaries 
in order to o btain the 2-D structure of an accretion disc (see 
lMarshll2005r ). 

The Doppler tomography method has been used by 
several authors f or the orbital solution in binary systems. 
Bagn uolo et al. I (1 19941 ') applied the Doppler tomography to 
a Be star (29 CMa in a binary system) for the first time. 
They used the tomography in order to produce separate UV 
spectra of the two stars and to obtain their spectral clas- 
sificat ions, mass ratio and so forth. Recently, IPeters et al. I 
(2008) applied the Doppler tomography method for recon- 
struction of UV spectra of the secondary FW CMa, because 
many authors had suggested that its variability in optical 
and UV spectra was due to the binary interaction although 
its binarity had been unknown. They also found that Ha 
and He I A6678 variations with its orbital period imply a 
hot region in the outer area of the Be disc of the primary, 
facing the secondary. However, there has never been a direct 
application of the tomography to the variability of the Be 
disc. 

We have used public software for the Doppler tomogra- 



phy method developed by Henk Sprufjf]. The V/R variation 
phase <&v/r (determined by the 500-day period) is adopted 
instead of the usually-assumed orbital phase, because it is 
not the binary motion but a one-armed density wave that ro- 
tates rigidly. Our data for the Doppler tomography consists 
of 47 Ha spectra (HJD 2453699.060 - HJD 2454902.965) 
and 37 U/3 spectra (HJD 2454412.184 - HJD 2454902.965). 
The coverage of the V/R variation is more than 2 cycles for 
the Ha lines, while < 1 c ycle for the H/3 l ines. The systemic 
velocity of —30 kms" - 1 jHutchingsl[l984l ) is taken into ac- 
count. However, the radial velocity of the Be star was not 
subtracted from the result, because it is very difficult to 
obtain it from each spectrum. Therefore, each line profile 
contains an uncertainty in wavelength. The absolute value 
of reporte d radial velocity ranges from several to several ten s 
of kms -1 l|Hutchings et al." I ll978l ; iGrundstrom et al. II2007T ) , 
so that each line profile contains an uncertainty of ~ 0.1 - 
1 A in wavelength. This does not affect the resulting map 
much (see below). 

Figure [9] shows results of the Doppler tomography. The 
results for the Ha line are displayed in the left column, and 
those for the H/? line in the right column. Top panel displays 
the obtained disc structures with the intensity in arbitrary 
unit in the velocity plane. At &v/r — 0, the observer is in 
the +a;-direction, and rotates clockwise. Obtained spectrum 
at that time is a projection of Doppler map along the V x 
axis. Note that the plus/minus sign of of V x and V y is de- 
fined in the x-y plane, and that if a structure, therefore, 
has a velocity in the — V^-direction, the observer obtains 
the spectrum with the redshift. The reconstructed spectra 
derived from the obtained disc structure and its deviations 
from the observed ones are shown in the middle and the bot- 
tom panels, respectively. The derived Doppler map shows a 
non-axisymmetric structure of the Be disc; a brighter re- 
gion is in V y > 0, while a fainter region in V y < 0. This 
feature is seen more clearly in the map of H/3 than in Ha. 
These maps suggest the presence of a global one-armed den- 
sity wave. Thus, our result provides a direct evidence of the 
one-armed oscillation in the Be disc. Interferometric obser- 
vations have alread y obtained such den sity enhancem ents in 
the br ight Be stars [Vakili et ail j 19981 ) for £ Tau and lBerid 
l|l999l ) for 7 Cas, bot h with GI2T; see also recent VLTI ob- 
servations of C Tau (IStefl et al. 1 120091 ) and its theoretical 
modeling (|Carciofi et al. I 20091 )]. However, this is the first 
time that a method other than the interferometry directly 
probed the two-dimensional distribution of the Be disc emis- 
sivity. 

There is 10 - 20 % difference between observed spectra 
and reconstructed ones. We should note here that 1) the 
radial velocity variations are not corrected in this analysis, 
as mentioned above, 2) the data are not sufficient especially 
in &v/r — 0.4 - 0.7, and 3) the non-axisymmetric structure 
due to the one-armed oscillation is possibly deformed even 
within one cycle. These may affect the resulting Doppler 
map to some extent. The first two points should be fixed 
by future observations, and more theoretical work should be 
done on the last point. 

In the case of isolated Be stars or binary Be stars with 
well-determined orbital elements, the structure revolving 



http://www.mpa-garching.mpg.de/~henk/ 
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Figure 8. Long-term variations of EW of the Ho (upper pan el) and H/3 (lower pa nel) lines : iHutchings et al~l Jl97ST ) (diam onds). lAab 
" 1985|1 (filled circles) iMotch et al. I il99lT) (fi lled triangles), IClark et aU jl998ah ( crosses) , Irlaigh et al. I j2004l l (vluses) , ICoe et al. 



3 (open ci^ l Gnmdstromrt trn a (sauares) , and this work (open triangles). 



with the period of V/R variation is seen more clearly than 
the current analysis. Hence, we can investigate whether the 
Doppler map contains a perturbation component. Subtract- 
ing the axisymmetric component from the Doppler map, we 
can obtain the perturbation component of the Be disc. Then, 
we can extend the investigation into the question that 1) the 
perturbation pattern is eccentric or spiral, and 2) one-armed 
oscillation is confined to inner part of the Be disc. If the os- 
cillation is confined to the inner part of the disc, the Doppler 
map in the velocity plane shows an axisymmetric structure 
only in the inner part, and the outer area of the map is 



non-axisymmetric. These studies will be achieved with the 
future work. 



As described above, the Doppler tomography method 
is very useful for examining the structure of the Be disc, 
although the information on time variability cannot be ob- 
tained unlike in interferometric observations. The current 
study showed that, for Be stars with the periodic V/R vari- 
ability the 2-D disc structure can be obtained by applying 
the Doppler tomography method to high-dispersion spectra 
with 2-m class telescopes. 
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Figure 9. Results from the Doppler tomography method applied to the Ha {left) and H/3 (right) line profiles. Top panel: Obtained disc 
structure in the velocity plane. The Unit of the structure is arbitrary. Middle panel: Reconstructed spectrum derived from obtained disc 
structure. Bottom panel: Deviation of the reconstructed spectrum from the observed spectrum. The Doppler tomography maps suggest 
that a non-axisymmetric structure exists in the Be disc of A0535+26. 



5 CONCLUSIONS 

We have carried out high-dispersion spectroscopic monitor- 
ing observations of the Be/X-ray binary A0535+26/V725 
Tau from November 2005 to March 2009. Regarding the 
long-term variability, our results are summarised as follows: 

(i) The orbital period is determined to be 110.2 days by 
the Fourier analysis of RXTE/ASM data. 

(ii) The EW(Ha) and EW(H^) indicate that the Be disc 
in A0535+26 has kept growing after the last giant outburst 
in 2005. 

(iii) From our observ ations, together with those by 
iGrundstrom et al. I (|2007l ). the period of the V/R variations 



is determin ed to be approxim ately 500 days, which is con- 
sistent with lReig et al. I (|2005h . 

(iv) The Doppler tomography method applied to the Ha 
and H/3 line profiles has revealed a non-axisymmetric disc 
structure precessing with the V/R period, which can be orig- 
inated from the global one-armed oscillation. 



Further observation is needed to improve these results. In 
particular, the observations of metallic absorption lines be- 
tween 4000 and 5000 A is very important in order to obtain 
the ephemeris and then the radial velocity of the Be star. 
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